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This  paper  is  a  discussion  on  the  estimation  of  impedance  parameters  of  H2/ air  fed  Proton  Exchange 
Membrane  Fuel  Cells  (PEMFC).  The  impedance  model  corresponds  to  the  Randles  electrical  equivalent 
circuit  accounting  for  charge  separation  and  transport  processes  in  the  cathode  catalyst  layer,  as  well  as 
for  oxygen  diffusion  through  the  backing  layer.  A  sensitivity  analysis  confirms  that  the  cathode  pa¬ 
rameters  are  not  correlated  and  that  the  consideration  of  the  anode  has  no  significant  impact  on  the 
estimation  of  their  values.  In  addition,  it  is  shown  that  the  diffusion  parameters  have  a  significant  impact 
in  the  low  frequency  domain  only,  at  least  with  this  model.  The  parameters  characterizing  charge  sep¬ 
aration  and  transport  processes  at  the  cathode  can  thus  be  estimated  with  the  high  frequency  impedance 
data,  independently  of  the  oxygen  transport  model.  Consequently,  even  in  the  absence  of  a  fully  validated 
oxygen  transport  impedance,  EIS  can  be  used  as  an  alternative  method  (to  classical  steady-state 
methods)  for  the  estimation  of  the  parameters  characterizing  the  cathode  reaction:  the  Tafel  slope  b, 
the  charge  transfer  coefficient  a  and  possibly,  the  exchange  current  density  jo.  This  reduces  significantly 
the  measuring  time  while  enhancing  the  accuracy  by  comparison  with  steady-state  methods. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Applied  to  fuel  cells,  Electrochemical  Impedance  Spectroscopy 
(EIS)  is  a  powerful  diagnosis  tool  that  can  be  used  during  operation 
with  minimal  perturbations  [1—3].  For  a  comprehensive  review 
about  the  application  of  EIS  to  Proton  Exchange  Membrane  Fuel 
Cells  (PEMFC)  we  refer  to  the  works  of  Yuan  et  al.  [4]  and  Niya  and 
Hoorfar  [5]  giving  an  overview  of  recent  advances  in  this  field. 
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EIS  allows  to  analyze  in  the  frequency  domain  the  different 
physical  and  chemical  processes  occurring  in  a  fuel  cell  by  imposing 
small  harmonic  perturbations  to  the  system  and  measuring  the 
impedance.  Consequently,  this  technique  can  be  used  for  in-situ 
and  ex-situ  studies  which  makes  it  a  polyvalent  and  accurate  tool 
by  comparison  with  most  of  the  classical  electrochemical  (steady- 
state)  methods.  A  common  application  of  EIS  is  the  identification  of 
the  origin  of  the  different  potential  losses  (i.e.,  ohmic,  kinetic  and 
mass  transport  losses)  and  their  estimation  [1,2, 6-9]  with  the  focus 
to  optimize  the  MEA  design  [10]  and  to  determine  optimal  working 
conditions  [1,11].  Beside  this,  it  is  possible  to  determine  the  elec¬ 
trical  properties  of  the  fuel  cell  components,  materials  and  in¬ 
terfaces  [12],  as  well  as  their  dependence  on  the  operating 
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conditions  [2,13,14]  and  on  aging  effects  like  carbon  corrosion  at  the 
cathode  [15  ,  catalyst  poisoning  by  carbon  monoxide  [14,16]  at  the 
anode  or  corrosion  of  the  flow  field  plates  [17].  EIS  can  also  be  used 
for  the  estimation  of  the  catalytic  active  area  [9],  although  other 
electrochemical  techniques  (e.g.,  cyclic  voltammetry)  are  more 
appropriate.  The  water  management  and  the  MEA  humidification 
playing  an  important  role  in  the  electrical  efficiency  of  a  fuel  cell, 
numerous  EIS  studies  focus  on  the  detection  of  dry  and  flooded 
conditions  and  on  their  influence  on  the  fuel  cell  operation  [1,9,14]. 

Obtaining  reliable  information  about  the  physical  processes  that 
occur  in  a  fuel  cell  is  not  a  simple  task  since  they  are  strongly 
coupled.  The  experimental  impedance  spectra  are  usually  inter¬ 
preted  in  terms  of  equivalent  electrical  circuits  (EEC)  consisting  of 
an  assembly  of  resistances  R,  capacitances  C  and  inductances  L. 
Almost  all  experimental  impedance  spectra  in  the  right  side  of  the 
complex  plane  can  be  fitted  with  an  EEC  provided  that  it  contains 
enough  R/C  elements  in  parallel.  As  a  consequence,  a  multitude  of 
more  or  less  sophisticated  EECs  can  be  found  in  the  literature  to 
describe  all  possible  forms  of  impedance  spectra.  However,  in  order 
to  identify  the  origin(s)  of  voltage  drop(s),  the  elements  of  the 
electrical  circuit  have  to  be  consistent  with  the  physical  phenom¬ 
ena  occurring  in  the  cell.  In  addition,  for  their  simultaneous 
independent  estimation,  the  impedance  parameters  have  to  be 
decorrelated.  Hence,  this  paper  does  not  aim  to  present  a  new 
impedance  model,  but  to  discuss  about  the  validity  of  the  estima¬ 
tion  of  the  parameters  that  characterize  charge  separation  and 
transfer  processes  in  H2/air  fed  PEMFC  when  using  EIS:  the  dis¬ 
cussion  is  based  on  numerical  and  experimental  studies  using  the 
simple  and  commonly  used  Randles  EEC  [18]. 

In  this  circuit,  charge  separation  and  transport  are  accounted  for 
by  a  high  frequency  resistance,  a  charge  transfer  resistance,  and  a 
double-layer  capacity.  The  oxygen  transport  through  the  porous 
media  at  the  cathode  is  modeled  by  a  finite  Warburg  element  [19], 
which  results  from  simple  assumptions:  one-dimensional  transport 
by  diffusion  only,  surface  reaction  and  a  constant  02  concentration  at 
the  GDL/air  channel  interface.  Since  the  oxygen  transport  parameters 
depend  strongly  on  the  chosen  impedance  model  [20-22],  we  focus 
on  the  estimation  of  those  characterizing  the  ORR.  It  is  shown  that 
their  values  can  be  estimated  independently  of  the  oxygen  transport 
impedance  through  the  analysis  of  the  impedance  spectra  in  the  high 
frequency  domain  only.  Some  results  of  the  estimation  of  parameters 
characterizing  the  cathode  reaction  are  presented,  i.e.,  the  Tafel  slope 
b  and  the  charge  transfer  coefficient  a  (or  the  exchange  current 
density  jo)  vs.  the  current  density.  Usually,  these  parameters  are 
determined  using  steady-state  methods  (e.g.,  the  Tafel  plot),  which 
leads  to  two  (or  three)  averaged  values  being  associated  with  the  low 
(the  intermediate)  and  the  high  current  density  regimes  [23-26]. 
Former  studies  show  that  EIS  can  be  used  as  an  alternative  to  the 
polarization  curve  for  the  determination  of  the  Tafel  slope  [24,25]. 
However,  in  these  studies,  the  values  of  b  are  averaged  over  the  low 
and  high  frequency  ranges.  In  contrast,  we  show  that  EIS  can  be  used 
to  estimate  not  only  b,  but  also  a,  and  possibly  jo  at  a  given  operating 
point  with  one  single  impedance  measurement.  In  addition,  their 
estimation  using  the  impedance  data  at  high  frequencies  only  is  in¬ 
dependent  of  oxygen  transport  losses,  which  characteristic  time- 
scales  correspond  to  much  lower  frequencies.  This  separation  of 
ORR  parameters  from  oxygen  transport  is  not  always  possible  when 
using  steady-state  methods.  EIS  could  thus  not  only  reduce  the 
measuring  time,  but  also  enhance  the  precision  of  the  estimation  of 
the  reaction  parameters  of  the  cathode. 

In  the  second  part  of  this  work,  we  present  results  of  a  sensitivity 
analysis  ensuring  the  absence  of  correlation  between  the  model  pa¬ 
rameters,  which  is  a  necessary  condition  for  their  reliable  estimation. 
More  generally,  the  sensitivity  analysis  allows  to  point  out  the  abilities 
and  limits  of  the  impedance  model  in  the  considered  operating 


conditions.  The  anode  is  usually  neglected  in  impedance  analyses 
because  the  voltage  drops  caused  by  this  electrode  are  significantly 
smaller  than  those  originating  from  the  cathode  [8,9,13,27,28].  This 
hypothesis  is  also  tested  through  the  sensitivity  analysis. 

2.  Theory 

The  expression  of  the  cell  impedance  that  is  used  for  the  pa¬ 
rameters  estimation  is  derived  in  the  following.  In  section  2.1,  the 
expressions  of  the  cell  potential  and  current  are  presented.  These 
expressions  are  used  in  Section  2.2  for  the  derivation  of  the  cell 
impedance. 

Fig.  1  shows  a  schematic  representation  of  the  MEA  geometry. 
The  structure  is  assumed  symmetric  to  the  membrane.  The  main 
hypotheses  of  the  model  are: 

•  The  fuel  cell  is  isothermal  and  isobaric. 

•  The  description  of  the  transport  phenomena  is  one-dimensional, 
the  x-axis  being  oriented  from  the  anode  to  the  cathode. 

•  The  diffusion  media  at  the  anode  and  the  cathode  are 
homogeneous. 

•  The  catalyst  layers  are  considered  as  additional  diffusion  media 
of  a  thickness  of  about  10  pm.  Consequently,  we  assume  a  sur¬ 
face  reaction  (at  x  =  0)  and  a  constant  oxygen  concentration  at 
the  GDL/air  channel  interface  (at  x  =  5). 

•  The  cell  being  fed  with  pure  hydrogen,  the  diffusion  limitations 
are  negligible  at  the  anode  side  and  cH2(x  =  0)  =  cf^. 

•  Oxygen  transport  is  assumed  to  occur  by  Fickian  diffusion  only. 

The  anode  is  taken  into  account  in  order  to  estimate  its  influence 
on  the  identification  of  the  cathode  impedance  parameters. 

Note  that  the  assumption  of  one-dimensional  mass  transport  is 
valid  for  small  area  MEAs  only.  In  a  more  general  case,  two 
dimensional  mass  transport  phenomena  like  oxygen  depletion  [20] 
and  concentration  oscillations  induced  by  the  measuring  signal 
propagating  along  the  cathode  gas  channel  [29-31]  have  to  be 
taken  into  account. 

2  A.  System  equations 

In  a  general  way,  the  cell  voltage  can  be  expressed  by 
^cell (0  =  Eo  ~  ^hf  xJcell(0  —  VcW  ~  Va (0?  O ) 

where  E0  denotes  the  open  circuit  voltage  (OCV)  and  Rhf  an  ohmic 
resistance  accounting  for  the  ionic  resistance  of  the  membrane  as 
well  as  for  the  electronic  resistances  of  the  various  MEA  compo¬ 
nents  (catalyst  layer,  GDL,  flow  field  plates).  rjc  and  13  a  stand  for  the 
voltage  losses  resulting  from  the  finite  rate  of  the  cathode  and 
anode  reactions,  respectively. 

With  reference  to  the  catalyst  area,  the  values  of  the  exchange 
current  density  at  the  cathode  j0)C  are  commonly  comprised  be¬ 
tween  10-6Acmpt2  and  10  9Acmpt2  [9,23,24,32].  Assuming  an 
electrode  roughness  y  =  100,  the  values  of  jo,c  vary  between 
10-4  A  cm~2cro  and  10~7  A  cm~2cro  with  reference  to  the  macro¬ 
scopic  area  of  the  MEA  (AmaCr0  =  Ap t/y),  which  is  significantly  lower 
than  the  current  densities  of  an  operating  cell.  In  addition,  the  02 
diffusivity  through  the  pores  of  the  GDL  and  catalyst  layer  is 
limited,  especially  when  the  cell  is  fed  with  air  instead  of  pure  02 
[28].  This  phenomenon  is  reinforced  by  a  possible  flooding  of  the 
gas  pores  in  the  presence  of  liquid  water.  The  oxygen  concentration 
is  consequently  lower  at  the  electrode  c02  (x  =  0,  t)  than  at  the 
GDL/air  channel  interface  c02(x  =  3,  t)  =  Cq2.  The  oxygen  reduc¬ 
tion  reaction  (ORR)  is  thus  irreversible  and  rjc  can  be  described  by  a 
Tafel  law  of  the  form 
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Fig.  1.  Schematic  representation  of  the  MEA  geometry. 


VdO  =  bc  In 


~  x  c52  ; 

jo.c  X  Cq2(X  =  0,t). 


which  implies  that  (/f  c)t  =  (/f  a>t  =  (/cell)t.  Note  that  the  current 
(2)  densities  refer  to  the  macroscopic  area  of  the  MEA. 


The  exchange  current  density  at  the  anode  jo,a  varies  generally 
between  5x1 0  4  A  cnip  2  and  3x1 0-2  A  cnip  2  in  acidic  media 
[24,33-35].  This  results  in  a  macroscopic  exchange  current  density 
comprised  between  0.05  A  cm~2cro  and  3  A  cm~2cro  (y  =  100), 
which  is  in  the  range  of  the  output  current  densities.  The  hydrogen 
oxidation  reaction  (HOR)  can  therefore  be  considered  as  reversible 
and  it  must  consequently  be  described  by  the  Butler-Volmer 
equation 


Jf,a(0  =Jo,a 


Ch2(X  =  Mexp^aW 


lh2 

cH+(x  =  0  ,t) 


;a,H; 


exp|  -5^) 


(3) 


2.2.  Cell  impedance 


In  the  dynamic  regime,  the  cell  impedance  is  defined  by  the 
ratio  of  the  variation  of  the  cell  potential  to  that  of  the  cell  current 
density  and  reads  as 


_  ^^cell 

-4/cell 


(7) 


In  order  to  derive  an  explicit  expression  of  Zce n,  the  variables 
have  to  be  decomposed  in  their  steady-state  and  dynamic  com¬ 
ponents.  Assuming  harmonic  oscillations  of  small  amplitude,  each 
variable  X  (X  =  ECeii  Jceii  Jf,m>  Vm,  co2)  can  be  decomposed  into  its 
steady-state  value  (X)t  and  a  sinusoidal  component  of  amplitude 
AX,  leading  to 


To  simplify  the  notations,  it  is  supposed  that  the  proton  trans¬ 
port  is  significantly  faster  than  the  reaction  rate,  which  implies  that 
cH+  (x  =  0,  t)  =  c^+.  When  working  with  pure  hydrogen,  mass 
transport  is  not  limiting  36]  and  we  consider  that  the  H2  concen¬ 
tration  at  the  electrode  is  equal  to  that  in  the  channel  c^2  •  Finally, 
we  assume  equal  values  for  the  Tafel  slope  of  the  two  half  reactions 
at  the  anode  (baHz  =  baH+  =  ba  [24]),  and  (3)  simplifies  to 


jf.a(t)  =  Jo, a 


exp 


w)- 


exp  - 


Va(0 


(4) 


The  parameters  bc,a  in  Eqs.  (2)  and  (4)  denote  the  Tafel  slopes  for 
the  two  half  reactions  at  the  cathode  and  anode  respectively  which 
are  defined  as 


be,  a 


RT 

ac,anF  ’ 


(5) 


where  ac>a ,  are  the  respective  charge  transfer  coefficients  and  n  the 
number  of  exchanged  electrons.  By  convention,  the  overpotential 
and  the  current  density  are  usually  counted  negative  at  the  cathode 
and  positive  at  the  anode.  However,  for  the  sake  of  simplicity,  we 
consider  only  positive  values  of  ijm,  jf,m  and  jo,m  at  both  electrodes 
(m  =  a,c). 

The  cell  current  density  jceii(t)  is  the  sum  of  the  faradaic  jf(t)  and 
the  capacitive  current  density  jc(t)  =  Cdi,m  x  m/0t 


JcellW  =Jf,m(t)  +  Cdl,m^,  (6) 

where  Cdi,m  stands  for  the  double-layer  capacity  at  electrode  m.  In 
steady-state,  the  capacitive  contribution  to  the  current  is  zero, 


Ecell(t)  =  (Ecell)t  +  A£cell  exp(lWt) 

Jcell(t)  =  {/cell > t  +  A/cell  exp(iwf) 

=  0'cell>t  +  %m  exp(iwt)  (8) 

Vm(t)  =  +  A?;m  exp(iwt) 

Co2(x,t)  =  (co2(x))t+  Ac02(x)exp(iwf). 

The  overbars  denote  complex  numbers.  The  amplitudes  of  the 
system  parameters  can  be  derived  from  (1)— (6)  which  results  in 


A^cell  —  ^hf  x  4/cell  ^Vc  ^ Va 

_  u  A/f,c  _  u  Acq2(0) 

—  bC(ycell)t  bC(co2(0))t  (9) 

A/f ,a  =  [exp  (^)  +  exp  (  -  ^)  ]  A?ja 
4/ cell  ~  4/f,m  +  i^Cd\,m  x  ^Vm- 

Using  (9)  in  (7)  allows  to  write  the  cell  impedance  in  an  explicit 
form 


-l 


Zee  11  =  Rhf  + 


1 


be _ be  a^q2  (0) 

y  (/cell )  t  (Co2(0))t  Ajfc 


Jo, a 


exp 


(Va)t 


exp  - 


+  iwCdic 


(Va)t 


iojC 


dl.a 


-1 


(10) 


To  simplify  the  notations,  we  introduce: 
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The  anode  charge  transfer  resistance 


i?cta  =  a  =  /a  _  Z±_ 


dJ  f,a  Aif.a 


exp 


<4a)t 


exp  - 


(rla)t 


The  cathode  charge  transfer  resistance 


(ii) 


The  Warburg  impedance  can  thus  also  be  expressed  as  a  function  of 
these  variables:  Zw(<5,  Deff).  The  effective  diffusivity  accounts  for  the 
geometrical  properties  of  the  medium  through  the  porosity  e  and 
the  tortuosity  t  and  for  the  presence  of  liquid  water  which  impacts 
their  values.  Deff  is  usually  expressed  in  the  form  of  Archie’s  law 
[38,39]. 

With  these  assumptions,  the  cell  impedance  Zce n  (10)  can  be 
written  as 


,  =  9^c  =  be 

bjf.c  (/cell)  t 

•  The  oxygen  transport  impedance 


(12) 


ZE  =  Rhf 


^w 


-1 


+  iwQi  c  J  +  (  — f  icoCc 


^dl,a 


-1 


(18) 


,  =  be  Acq2(0) 

(cO2(0))t  Ajf 


(13) 


An  appropriate  oxygen  transport  model  is  required  to  obtain  an 
explicit  expression  of  Zc 02,  because  the  estimation  of  the  oxygen 
diffusion  parameters  (Deff,<5)  depends  strongly  on  the  transport 
phenomena  taken  into  account  by  the  model.  A  detailed  study 
about  oxygen  transport  impedance  expressions  would  be  beyond 
the  scope  of  this  work.  For  more  information,  see  for  example  [20- 
22,29,30].  In  the  present  study,  oxygen  transport  is  accounted  for 
by  a  finite  Warburg  impedance  Zw  allowing  a  first  rough  estimation 
of  Deff  and  5  [37].  The  Warburg  impedance  (14)  is  obtained  by 
solving  the  Fick  diffusion  equations  in  the  dynamic  regime  with  the 
boundary  conditions  depicted  in  Fig.  1  [19].  This  corresponds  to  a 
purely  diffusive  one-directional  oxygen  transport  through  homo¬ 
geneous  porous  media  and  Zw  is  given  by 


Zw  =  R{ 


tanh^iwid 

1 


(14) 


where  the  diffusion  resistance  R d  corresponds  to  the  low  frequency 
limit  of  Zw 


Rd  =  lim  =  - - 7-, 

d  4FDeff(co2(0))t 


(15) 


and  id  is  the  characteristic  diffusion  time 


(16) 


The  time-averaged  oxygen  concentration  at  the  reaction  inter¬ 
face  (cO2(0))t  is  obtained  by  solving  the  Fick  equation  in  steady- 
state,  leading  to  [20] 


<co2(0)>t  =  c52  —  (17) 

It  is  worth  mentioning  that  Rd  and  rd  depend  both  on  the 
diffusion  thickness  <5  and  effective  oxygen  diffusion  coefficient  Deff. 


Fig.  2  represents  the  cell  impedance  described  by  Eq.  (18)  in 
terms  of  equivalent  electrical  circuits.  The  anode  EEC  is  composed 
by  a  charge  transfer  resistance  Rc t>a  in  parallel  with  a  double-layer 
capacity  Cdi>a  (dashed  lines).  The  membrane  and  the  cathode  are 
represented  by  Randles  EEC  [18  including  a  high  frequency  resis¬ 
tance  Rhf.  a  charge  transfer  resistance  Rct,c»  a  double-layer  capacity 
Cdi)C  and  a  finite  Warburg  impedance  Zw  [19]  (solid  lines). 


3.  Experimental  setup 

The  experimental  results  are  obtained  with  the  mono-channel 
single  cell  presented  in  Fig.  3.  The  serpentine  flow  fields  (of 
0.7  x  1.0  mm2  in  section  and  33  cm  in  length)  are  symmetric  at  the 
anode  and  cathode.  Although  only  global  impedance  data  are 
treated  in  this  work,  it  must  be  noted  that  the  cell  was  designed  to 
perform  also  locally  resolved  EIS  [37]. 

The  MEA  (Johnson  Matthey)  with  a  flat  area  Amacro  =  7.87  cm2 
consists  of  a  PFSA  (Gore)  polymer  membrane  of  thickness 
<5m  =  30  |im  and  active  layers  (<5ai  ~  10  fim)  with  an  average  Pt 
loading  of  0.4  mgPt  cm-2.  The  GDL  (Toray™  TGP-EI-060)  is  190  ftm 
thick  with  a  porosity  (without  compression)  e  =  0.78.  The  fuel  cell  is 
fed  in  counter-flow  (at  1  atm)  with  dry  hydrogen  and  humidified  air 
with  the  corresponding  stoichiometry  factors:  Sh2  =1.2  and 
Sair  =  3.  The  air  is  humidified  ( RH  ~  78%)  and  pre-heated  thanks  to 
a  bubbler.  The  connecting  pipes  are  slightly  heated  in  order  to  avoid 
condensation  between  the  outlet  of  the  bubbler  and  the  air  channel 
inlet.  The  cell  temperature  is  controlled  (Tceii  =  56  ±  1°C)  by  an 
electric  heating  system  fixed  to  the  anode  bipolar  plate.  Thermal 
gradients  exist  probably  within  the  cell  because  of  the  thermal 
resistances  of  the  various  components  and  those  of  the  interfaces 
between  them  [40]. 

The  impedance  measurements  are  performed  in  galvanostatic 
mode  at  frequencies  ranging  logarithmically  from  0.025  EIz  to 
500  EIz  with  a  resolution  of  10  points  per  decade.  The  peak-to-peak 
sinusoidal  perturbation  is  comprised  between  5%  and  10%  of  the 
cell  current  (at  (/cell)t  =  0.5  A  cm-2,  25  mA  cm-2  <  Ajcell  < 
50  mA  cm-2  and  31  mV  <  AEcell  <  62  mV).  It  was  checked  that  in 
this  range,  the  AC  signal  amplitude  does  not  affect  the  impedance 
measurements. 


Fig.  2.  MEA  equivalent  electrical  circuit(s).  Dashed  lines:  anode  equivalent  circuit  used  for  sensitivity  analysis  only.  Solid  lines:  Randles  equivalent  circuit  [18]  of  the  cathode  and 
electrolyte  used  for  parameters  estimation. 
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Fig.  3.  Mono-channel  single  cell  with  a  flat  area  Amacro  —  7.87  cm2.  The  current  col¬ 
lectors  form  the  serpentine  flow  field. 

The  data  acquisition  and  processing  system  is  depicted  in  Fig.  4. 
Since  the  cell  impedance  is  almost  entirely  determined  by  the 
cathode,  the  membrane  and  some  contact  resistances,  it  is  possible 
to  use  the  anode  as  a  reference  electrode  [9].  A  passive  electronic 
circuit  is  used  to  measure  the  current  intensities  with  the  least 
possible  perturbations:  each  of  the  current  collectors  is  connected 
to  a  10  mQ  shunt  resistance.  The  whole  experimental  setup  is 
controlled  by  a  Labview  program. 

The  impedance  parameters  are  determined  by  fitting  the 
experimental  spectra  with  Eq.  (18)  using  a  complex  nonlinear  least 
square  method.  In  practice,  we  use  the  Matlab  function  fminsearch 
that  determines  the  minimum  of  the  function  AZ,  defined  as  the 
square  of  the  difference  between  the  experimental  Zexp  and  the 
theoretical  impedance  Zceii  (18)  —  A Z=  ( Zexp  -  Zceii)2,  starting  from 
an  initial  set  of  model  parameters  (Rhf,  Rct,c,  Qi.o  Rd,  t,  Rct,a,  Cdi.a)- 
fminsearch  uses  the  Nelder-Mead  simplex  algorithm  of  Lagarias 
et  al.  [41]. 

4.  Results  and  discussion 

In  Section  4.1,  we  present  and  discuss  some  results  of  the  pa¬ 
rameters  estimation  characterizing  the  charge  separation  and 
transfer  processes  at  the  cathode  obtained  using  the  impedance 
model  of  the  Randles  ECC,  i.e.,  the  Tafel  slope  b,  the  charge  transfer 
coefficient  a  and  the  exchange  current  density  jo.  These  parameters 
can  be  determined  with  the  charge  transfer  resistance  Rct 


independently  of  the  oxygen  transport  model  by  fitting  the  spectra 
in  the  high  frequency  domain  only  (cf.  sensitivity  analysis  in  Section 
4.2).  It  is  worth  mentioning  that  the  values  of  b  obtained  with  Rct 
depend  on  the  reaction  kinetics  and  on  the  proton  conduction 
through  the  electrode  volume  ( lumped  b ).  As  the  cathode  param¬ 
eters  only  are  considered,  the  subscripts  relative  to  the  electrodes  (a 
or  c)  are  omitted  in  this  part  to  simplify  the  notations. 

A  complementary  sensitivity  analysis  relative  to  the  parameters 
of  the  impedance  expression  in  Eq.  (18)  is  carried  out  in  Section  4.2 
in  order  to  investigate  the  domain  of  the  validity  of  the  Randles  EEC 
(Fig.  2)  used  for  the  estimation  of  the  cathode  parameters. 

4.1.  Estimation  of  cathode  reaction  parameters 

The  following  experimental  data  correspond  to  mean  values 
resulting  from  a  series  of  five  measurements  carried  out  under 
identical  conditions.  The  corresponding  standard  deviations  are 
also  shown  in  the  figures.  The  cell  was  operated  as  described  in 
section  3,  at  current  densities  jceii  comprised  between  0.013  A  cm-2 
and  1.1  A  cm-2. 

Fig.  5  depicts  the  charge  transfer  resistance  Rct  as  a  function  of 
jceii-  The  small  error  bars  indicate  the  good  reproducibility  of  the 
estimation  using  the  Randles  EEC.  The  uncertainties  increase  at  low 
current  densities.  This  is  related  to  the  higher  amplitude  of  the  AC 
signal  (with  reference  to  the  DC  signal)  required  to  perform  the 
impedance  measurements  at  low  current  densities. 

Rct  remains  uniform  over  the  major  part  of  the  current  density 
range.  However,  it  increases  steeply  at  low  current  densities 
(jceii  <  0.1  A  cm-2),  which  is  consistent  with  the  steep  increase  of 
the  activation  overpotential  in  this  range.  At  very  low  current 
densities  (near  the  OCV),  the  ORR  can  be  described  by  a  linearized 
form  of  the  Butler— Volmer  law  which  results  in  combination  with 
the  definition  of  Rct  at  the  cathode  (12)  in 


A  more  detailed  derivation  of  the  above  expression  is  given  in 
the  Appendix.  It  must  be  noted  that  Rct  is  related  to  the  global  ORR 
per  mole  of  produced  water  (l/202+2H++2e-  -►  H2O)  so  that  n  is 
set  to  2  in  our  case.  Nevertheless,  the  expressions  of  the  Tafel  slope 
and  of  the  charge  transfer  resistance  derived  above  remain  valid  no 
matter  the  number  of  exchanged  electrons. 

The  values  of  jo  are  generally  comprised  between  10  4  A  cm-2 
and  10-7  A  cm-2  [9,23,24,32]  with  reference  to  the  MEA  flat  surface 
(with  y  =  100).  Starting  from  these  results,  the  value  of  the  charge 
transfer  resistance  near  the  OCV  could  be  (19): 


Fig.  4.  Data  acquisition  and  processing  setuP  for  the  imPedance  measurements. 
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Fig.  5.  Cathode  charge  transfer  resistance  Rct  identified  from  El  spectra  using  the 
Randies  equivalent  circuit  vs.  cell  current  density  jceii. 


j0  =  10~4  A  cm~2  — >Rct(OCV)  =  71  Q  cm2 
jo  =  10“7  Acm-2^Rct(OCV)  =  7.1  104ficm2. 

The  behavior  of  Rct  obtained  with  our  EIS  measurements  at  low 
current  densities  is  thus  consistent  with  the  expectations, 
notwithstanding  the  significant  variations  of  the  data  available  in 
the  literature.  However,  the  estimation  of  the  electrode  reaction 
parameters  ( b  and  ac)  must  be  made  at  current  densities  where  a 
Tafel  law  (2)  can  be  used.  Consequently,  the  values  of  Rc t  estimated 
at  the  lowest  experimental  current  densities  (jceii  <  0.013  A  cnrT2) 
are  not  considered  in  the  following. 

The  most  commonly  used  method  for  the  estimation  of  the  Tafel 
slope  is  the  Tafel  plot  (plot  of  the  cathode  overpotential  vs.  the 
common  logarithm  of  the  current  density)  yielding  a  straight  line 
of  slope  b  10  (expressed  in  V  dec-1)  and  a  y-intercept  at 
-b10  log(/0Co2/Co2).  This  method  leads  to  two  (or  sometimes  three) 
average  values  of  bio1  for  the  oxygen  reduction  reaction  (in  Pt/C 
electrodes):  at  low  current  densities  (E  >  0.8  V),  bio  ~  60  mV  dec  1 
and  at  high  current  densities  (E  <  0.8  V),  bio  ~  120  mV  dec-1 
[8,23,25,42-44].  The  change  in  the  Tafel  slope  is  generally 
explained  by  a  shift  in  the  mechanism  from  a  reduction  on  a  Pt- 
oxide  covered  surface  (at  low  current  densities)  to  a  reaction  on 
an  oxide  free  surface  (at  high  current  densities)  characterized  by 
different  adsorption  isotherms  (Langmuir,  Temkin)  [25,44].  It 
should  be  noted  that  these  values  are  only  orders  of  magnitude  and 
the  results  reported  in  the  literature  vary  usually  in  a  range  of  about 
20  mV  dec-1  depending  on  temperature,  gas  pressure,  oxidant  gas 
(pure  O2  or  air)  and  investigated  system,  i.e.,  single  electrode,  single 
cell  or  stack  [9,23-25,42,43].  For  instance,  by  fitting  the  Tafel  plots 
of  a  PEMFC  operated  at  25  °C  with  H2  and  air  at  P  =  1  atm,  Ciureanu 
et  al.  9]  measured  a  Tafel  slope  varying  from  about  130  mV  dec-1  to 
270  mV  dec-1.  Regardless  of  the  exact  values  of  bio,  all  studies  in  the 
literature  report  consistently  an  increase  (doubling)  of  the  Tafel 
slope  from  the  low  to  the  high  current  density  regime. 

Using  EIS,  the  Tafel  slope  (relative  to  the  MEA  macroscopic  area) 
can  be  determined  directly  along  with  the  charge  transfer  resis¬ 
tance:  b  =  Rct  x  jceii  (12).  This  presents  the  advantage  of  identifying 
the  Tafel  slope  at  a  given  operating  point  by  one  single  El 
measurement. 

Fig.  6  depicts  the  Tafel  slope  (in  mV  dec-1)  derived  from  the 
charge  transfer  resistance  (Fig.  5)  as  a  function  of  jCeii*  bio  increases 
continuously  from  23  mV  dec-1  at  jceii  =  0.063  A  cm-2  to 


1  With  b10  =  In  10  x  b  =  2.3  x  b. 


384  mV  dec  1  at  jceii  =  1.08  A  cm-2.  The  resulting  values,  as  well  as 
the  increase  with  the  current  density  are  in  good  agreement  with 
former  studies  [8,23,25,42-44].  However,  instead  of  the  two  slopes 
regime  that  is  usually  reported,  bio  increases  continuously.  For 
comparison,  we  also  determined  the  value  of  bio  from  the  Tafel  plot 
and  obtained  bi0  =  67  mV  dec-1.  This  leads  to  the  conclusion  that 
the  Tafel  plot  yields  a  mean  value  of  bi0  averaged  over  the  experi¬ 
mental  current  density  range. 

The  Tafel  slope  can  related  to  the  (cathode)  charge  transfer  co¬ 
efficient  ac  using  Eq.  (5).  Considering  the  global  ORR  reaction  where 
n  =  2  this  yields  b  =  RT/(ac  x  2 F).  The  charge  transfer  coefficient  is 
an  important  parameter  characterizing  the  electrochemical  activity 
of  an  electrode.  The  usual  method  consists  in  estimating  this 
parameter  from  the  polarization  curve.  The  results  reported  in  the 
literature  (cf.  Table  1)  show  a  large  dispersion  of  the  values 
[25,26,45,46],  all  the  more  so  since  ac  depends  on  the  temperature 
in  the  high  slope  regime  25]. 

In  the  majority  of  the  cases,  this  parameter  is  comprised  be¬ 
tween  1  and  0.4  and  decreases  with  the  current  density.  For 
modeling  purpose,  ac  is  often  set  to  0.5,  which  corresponds  to 
symmetrical  forward  and  backward  reactions  at  the  cathode. 
However,  its  actual  value  depends  strongly  on  the  operating  con¬ 
ditions.  Measurements  performed  by  Parthasarathy  et  al.  [25]  show 
a  decrease  of  ac  from  1  (low  Tafel  slope  regime)  to  0.5  (high  Tafel 
slope  regime).  Danilov  et  al.  [26]  identify  3  different  values  of  the 
charge  transfer  coefficient  for  low,  intermediate  and  high  current 
densities:  their  results  confirm  the  decrease  of  ac  with  increasing 
current  density. 

Fig.  7  depicts  the  values  of  ac  deduced  from  the  charge  transfer 
resistance  in  Fig.  5  as  functions  of  the  current  density  jceii-  The 
uncertainty  on  the  estimation  of  ac  increases  at  low  current  den¬ 
sities,  because  of  the  lower  accuracy  of  the  impedance  measure¬ 
ments.  However,  regardless  of  the  lower  precision,  ac  approaches  1 
at  low  current  densities  and  decreases  with  increasing  cell  current. 
The  charge  transfer  coefficients  determined  with  Rct  as  well  as  their 
variation  with  the  current  density  are  in  good  agreement  with 
those  observed  in  former  studies  (cf.  Table  1)  confirming  that  EIS 
can  be  used  as  an  alternative  to  Tafel  plot.  Actually,  EIS  could 
possibly  provide  more  precise  estimates  than  the  steady-state 
methods,  because  ac  can  be  determined  directly  for  a  given  oper¬ 
ating  point. 

As  mentioned  above,  the  different  forms  of  oxygen  adsorption 
corresponding  to  different  adsorption  isotherms  influence  the 
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Fig.  6.  Tafel  slope  bw  derived  from  the  charge  transfer  resistance  in  Fig.  5  vs.  jceii- 
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Table  1 

Cathode  charge  transfer  coefficient  ac  as  given  in  the  literature. 


T  1°C] 

Slope  regime 

ac 

Ref. 

70 

Low  slope 

0.95 

[25] 

High  slope 

0.60 

70 

Low  slope 

0.89 

[26] 

High  slope 

0.43 

80 

Low  slope 

1 

[32] 

High  slope 

0.61 

reaction  mechanism  [25,44]  which  is  lumped  in  the  Tafel  slope.  For 
a  more  precise  estimation  of  ac  the  change  in  the  reaction  mech¬ 
anism  should  be  taken  into  account  through  the  value  of  n 
according  to  the  predominant  reaction  mechanism.  Furthermore, 
note  that  the  charge  transfer  resistance  and  the  diffusion  imped¬ 
ance  are  correlated  at  intermediate  frequencies,  which  is  shown  by 
the  sensitivity  analysis  presented  in  the  following  section.  This 
correlation  can  also  be  put  forward  using  oxygen  transport  models 
alternative  to  the  classical  Warburg  approach:  the  values  of  Rct  (and 
thus  of  bio  and  ac)  change  slightly  when  the  parameters  estimation 
is  performed  with  more  sophisticated  one-  and  pseudo  two- 
dimensional  oxygen  transport  impedances  (accounting  for 
convective  effects  along  and  perpendicular  to  the  electrode,  as  well 
as  for  a  finite  air  stoichiometry  [20—22]).  However,  this  does  not 
affect  the  behavior  of  the  electrode  reaction  parameters  as  a  func¬ 
tion  of  the  current  density  and  the  results  presented  above  remain 
qualitatively  valid. 


4.2.  Sensitivity  analysis 

The  following  sensitivity  analysis  aims  to  verify  the  reliability  of 
the  ECC  of  Fig.  2  for  the  estimation  of  the  cathode  impedance  pa¬ 
rameters.  To  this  end,  the  following  issues  are  considered: 


All  sensitivity  coefficients  Cx  are  normalized  according  to 
X  dZcell 

~  llZcdlll  dX  ’ 


(20) 


with  X  =  Rhf ,  Rc tiD  Cdi.o  Rd,  Td.  Rc t,a,  Qi.a  and  Zceii  standing  for  the 
global  impedance  of  the  MEA,  as  defined  by  Eq.  (18).  Cx  is  an  in¬ 
dicator  of  the  sensitivity  of  the  cell  impedance  to  the  parameter  X 
with  absolute  values  varying  between  0  and  1.  It  is  generally 
considered  that  when  the  absolute  value  of  Cx  is  lower  than  0.05, 
the  impact  of  X  on  the  cell  impedance  is  lower  than  the  experi¬ 
mental  data  scattering.  In  this  case,  a  consistent  parameter  esti¬ 
mation  is  not  possible  anymore. 

The  sensitivity  coefficients  Cx  are  complex  functions  of  the 
model  parameters.  In  what  follows,  the  Cx  are  obtained  using  the 
impedance  parameters  that  were  estimated  in  the  operating  con¬ 
ditions  given  in  section  3  at  jceii  =  0.5  A  cm-2  (Table  2). 


4.2 A.  Sensitivity  analysis  with  reference  to  the  cathode  parameters 

Fig.  8(a)  and  (b)  show  the  real  and  imaginary  parts  of  the 
sensitivity  coefficients  Cx  relative  to  the  cathode  parameters  as 
functions  of  the  frequency  v  in  a  range  comprised  between  1  mHz 
and  10  kHz.  The  absolute  values  of  Cx  may  vary  in  other  operating 
conditions.  However,  the  curves  in  Fig.  8  reflect  well  the  relative 
influence  of  the  different  model  parameters  over  the  whole  fuel  cell 
operating  domain. 

It  can  be  seen  that  the  sensitivity  coefficients  are  significant  in 
the  frequency  range  0.1  Hz  -  1  kHz.  For  frequencies  higher  or 
lower  than  this  interval,  the  values  of  Cx  are  either  negligibly 
small  (Cx  <  0.05)  or  constant  for  several  parameters  indicating 
that  correlations  may  exist  between  them.  Measurements  per¬ 
formed  with  frequencies  beyond  this  interval  would  thus  not 
contribute  to  a  better  precision.  The  observations  relative  to  the 
sensitivity  coefficients  of  the  cathode  parameters  can  be  sum¬ 
marized  as  follows: 


•  Possible  correlations  between  the  model  parameters. 

•  Possibility  to  estimate  the  cathode  parameters  without 
modeling  the  anode. 

•  Impact  of  the  low  frequency  parameters  (oxygen  transport 
impedance)  on  the  estimation  of  the  intermediate  and  high 
frequency  parameters  (fthf,  Rct  and  Cdi). 


Fig.  7.  Cathode  charge  transfer  coefficient  ac  deduced  from  EIS  measurements  vs.  jceii- 


•  The  sensitivity  coefficient  of  the  high  frequency  resistance  CRhf 
(black  solid  line  in  Fig.  8(a))  is  a  real  number.  R^  is  the  parameter 
with  the  highest  influence  on  the  cell  impedance  for  frequencies 
above  10  Hz:  its  sensitivity  coefficient  is  higher  than  0.4. 

•  The  real  and  imaginary  parts  of  the  sensitivity  coefficients 
relative  to  Rc t)C  (green  dotted  lines)  and  Cdi,c  (purple  dashed 
lines)  reach  their  maximum  (absolute)  values  between  10  Hz 
and  1  kHz.  In  this  range,  CRCt,c  and  Ccdi.c  climb  to  about  0.4. 

•  The  sensitivities  of  the  diffusion  impedance  parameters  -  CRd 
(pink  dash-dot  lines)  and  CT d  (blue  solid  lines)  -  reach  their 
highest  values  at  frequencies  below  10  Hz.  These  values  are 
generally  lower  than  those  of  the  cathode  reaction  parame¬ 
ters:  the  maximum  of  CRd  is  about  0.3  and  CXd  does  not  exceed 
0.1. 

These  results  confirm  that  the  time-scales  of  mass  transport 
and  of  charge  separation  and  transfer  are  clearly  separated  in  this 
impedance  model.  R hf,  Rct,c  and  Cdi,c  can  therefore  be  called  high 
frequency  parameters,  whereas  the  parameters  of  the  diffusion 
impedance  (d  and  Rd)  are  referred  to  as  the  low  frequency  pa¬ 
rameters.  These  results  also  put  forward  that  R^,  Rc t)C  and  Cdi,c  can 
be  estimated  without  the  use  of  an  oxygen  transport  impedance 
when  fitting  the  high  frequency  part  of  the  spectra  only.  Note 
that  this  result  is  valid  for  the  EEC  in  Fig.  2  only:  sensitivity  an¬ 
alyses  are  specific  to  each  model  and  should  be  carried  out 
systematically. 

It  has  to  be  noted  that  the  sensitivity  coefficients  relative  to  Rc t)C 
and  Rd  show  a  similar  behavior  at  intermediate  frequencies  indi¬ 
cating  that  these  parameters  could  be  correlated  over  this 
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Table  2 

Mean  cathode  and  anode  impedance  parameters  with  their  corresponding  standard  deviations  estimated  from  a  series  of  5  impedance  spectra  (Fig.  10)  and  anode  parameters 
as  they  can  be  found  in  the  literature  47  (third  row). 


Rhf  [fi  cm2] 

Rct, c  cm2] 

Cdi,c  [mF  cm  2] 

Rd  [Q  cm2] 

Td[s] 

Ret, a  [A  cm2] 

Cdi,a  [F  cm  2] 

0.156  ±  0.2% 

0.153  db  2% 

0.148  ±  0.06% 

0.27  ±  0.6% 

0.23  ±  0.8% 

0.23  db  0.8% 

22.4  ±  0.1% 

19.7  ±  0.5% 

19.7  db  0.1% 

0.158  ±  2% 
0.157  db  1% 

0.066  db  3% 

0.066  ±  3% 

0.005  db  20% 
0.007a 

0.2  db  50% 

0.75a 

1 

0.97a 

a  Fit  results  obtained  by  Meland  et  al.  [47]  on  a  cell  fed  with  hydrogen  at  both  electrodes  at  T  =  50  °C  and  P  =  1  atm. 


frequency  range.  It  is  thus  possible  that  the  choice  of  the  oxygen 
transport  model  affects  the  estimation  of  the  charge  transfer 
resistance  when  using  the  Randles  EEC  The  impact  of  this  possible 
correlation  on  the  estimation  of  the  high  frequency  parameters  is 
treated  with  more  details  below. 


4.2.2.  Sensitivity  analysis  with  reference  to  the  anode  parameters 

For  the  sake  of  clarity,  the  anode  sensitivity  coefficients  are 
plotted  separately  in  Fig.  9  as  functions  of  the  frequency.  Similarly 
to  the  cathode  sensitivity  parameters,  Ccdi.a  (red  and  pink  solid 
lines)  and  CRCt,a  (dark  and  light  green  dashed  lines)  reach  their 
maximum  (absolute)  values  at  frequencies  between  10  Hz  and 
1  kHz,  but  at  slightly  higher  frequencies  than  the  corresponding 
cathode  sensitivities.  However,  the  absolute  values  of  the  anode 
sensitivities  do  not  exceed  0.01,  whereas  the  cathode  coefficients 
reach  about  0.4.  The  sensitivity  of  the  cell  impedance  relative  to  the 
cathode  parameters  is  about  40  times  higher  than  to  those  related 
to  the  anode.  This  leads  to  the  conclusion  that  the  anode  can  be 
neglected  in  the  impedance  model,  which  is  consistent  with  results 
from  previous  articles  [8,9,13,27,28].  It  has  to  be  noted  that  this 
conclusion  is  valid  only  when  the  anode  operates  normally  (i.e.  in 
the  absence  of  phenomena  such  as  fuel  starvation,  water  flooding 
or  electrode  degradation). 

As  fuel  cell  impedance  spectra  are  often  depressed  semicircles, 
Constant  Phase  Elements  (CPE)  are  sometimes  used  to  enhance  the 
quality  of  the  fitting  curves  [9,47-49].  There  exist  various  expres¬ 
sions  of  the  impedance  of  a  constant  phase  element  Zcpe  [48]. 
However,  in  most  of  the  cases,  ZCpe  is  written  as 


Zcpe 


1 


Ql,a(IW) 


a  j 


(21) 


where  a  <  1.  This  impedance  is  purely  resistive  for  a  =  0,  capacitive 
for  a  =  1  and  inductive  for  a  =  -1;  generally  a  e  [0.5,1].  The 
physical  meaning  of  CPE  is  usually  related  to  the  distribution  of  the 
capacitance  through  the  catalyst  layer  [9,50,51].  The  impact  of  the 
CPE  in  the  anode  equivalent  circuit  is  analyzed  in  Fig.  9:  Ca  (light 
and  dark  blue  dash-dot  lines)  not  only  governs  the  anode  imped¬ 
ance  over  the  whole  frequency  range,  but  it  has  a  similar  behavior 
as  Ccdi.a*  indicating  a  possible  correlation  between  both  parameters. 
A  simultaneous  estimation  of  the  anode  parameters  is  (most 
probably)  not  possible  when  using  a  CPE  instead  of  a  capacitance  in 
this  EEC.  This  could  be  confirmed  experimentally:  the  fitting  curves 
do  not  converge  when  using  a  CPE  in  the  anode  EEC:  the  exponent  a 
(at  least)  has  to  be  imposed  as  a  fixed  model  parameter.  In  order  to 
avoid  problems  related  to  the  use  of  a  CPE,  a  simple  capacitance  was 
used  in  the  anode  EEC. 


4.2.3.  Experimental  validation  of  the  numerical  sensitivity  analysis 
The  results  of  the  (numerical)  sensitivity  analysis  can  be 
confirmed  by  fitting  experimental  spectra  with  the  equivalent 
circuit  of  Fig.  2.  Fig.  10  shows  a  mean  El  spectrum  (average  values 
and  standard  deviations)  resulting  from  a  series  of  5  spectra 
measured  under  identical  conditions  (those  of  the  sensitivity 
analysis,  cf.  Section  3).  For  frequencies  below  0.4  Hz,  the  data 
exhibit  an  important  scattering  which  can  be  linked  to  the 
evacuation  of  water  droplets  through  the  air  channel.  Impedance 


Fig.  8.  Real  (a)  and  imaginary  part  (b)  of  the  cathode  sensitivity  coefficients  Cx  as  functions  of  the  frequency  v  for  the  set  of  impedance  parameters  in  Table  2;  X  =  Rhf,  Rct,c.  Qi,c>  ^d. 

^d>  ^ct,a>  Oil, a- 
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Fig.  9.  Real  and  imaginary  parts  of  the  anode  sensitivity  coefficients  as  functions  of  the  frequency  using  the  set  of  impedance  parameters  in  Table  2. 


data  measured  at  frequencies  below  0.4  Hz  are  thus  not  retained 
for  the  parameter  estimation.  Furthermore,  this  avoids  compli¬ 
cations  resulting  from  possible  inductive  features  sometimes 
observed  at  low  frequencies  [8,36,52-54].  The  experimental 
spectra  exhibit  also  sometimes  an  inductive  loop  in  the  kHz 
regime,  which  is  consistent  with  previous  observations  [9,54,55]. 
Our  analyses  show  that  the  sensitivity  coefficients  reach  signifi¬ 
cant  values  between  0.4  Hz  and  200  Hz  (Fig.  8)  and  that  in 
addition,  the  model  parameters  do  not  seem  to  be  correlated  in 
this  range.  Consequently,  parameter  estimations  are  limited  to 
this  interval  where  inductive  loops  are  without  remarkable 
consequences.  This  allows  also  to  exclude  the  straight  line  at  45° 
that  appears  in  the  kHz  regime  and  which  is  linked  to  protonic 
resistance  of  the  volumetric  cathode  [1,56,57]. 

Table  2  shows  the  estimated  mean  cathode  and  anode  imped¬ 
ance  parameters  (leading  to  the  theoretical  impedance  spectrum  in 
Fig.  10;  red  curve)  with  their  expected  standard  deviations.  In  the 
first  row,  we  present  the  cathode  parameters  obtained  using  the 
Randles  EEC  only.  They  are  consistent  with  values  commonly  found 
in  the  literature  [27,58,59].  In  addition,  their  corresponding  stan¬ 
dard  deviations  do  not  exceed  3%,  which  is  an  indicator  of  the 
validity  of  this  parameters  estimation. 

The  second  row  of  Table  2  shows  the  parameters  resulting  from 
the  estimation  using  the  anode  and  cathode  EEC  (Fig.  2).  Actually, 
the  direct  estimation  of  the  8  parameters  is  not  possible  because  of 
cross-correlations  between  them.  Therefore,  the  following  method 
was  employed: 


Fig.  10.  Experimental  mean  spectrum  resulting  from  a  series  of  5  measurements  under 
identical  conditions  and  theoretical  spectrum  computed  after  identification  of  the 
model  parameters. 


1.  Estimation  of  the  cathode  parameters  keeping  constant  the 
anode  parameters. 

2.  Estimation  of  the  anode  parameters  keeping  constant  the 
former  estimated  cathode  parameters. 

3.  Re-estimation  of  the  cathode  parameters  keeping  constant 
again  the  anode  parameters. 

The  values  of  the  cathode  parameters  exhibit  almost  no  varia¬ 
tion  between  steps  1  and  3.  It  can  thus  be  confirmed  that  the 
estimation  of  the  anode  parameters  has  no  significant  influence  on 
the  cathode  parameters.  Actually,  the  only  difference  lies  in  the 
value  of  Rhf-  when  considering  the  cathode  only,  Rhf  corresponds  to 
the  sum  of  the  high  frequency  resistance  and  the  anode  charge 
transfer  resistance  Rcta  resulting  from  the  estimation  when 
considering  both  electrodes. 

Note  that  although  our  experimental  values  of  Rct,a  and  Cdi,a 
exhibit  a  scattering  reaching  50%,  they  are  in  good  agreement  with 
those  encountered  in  the  literature  for  similar  operating  conditions 
[47  (third  row  of  Table  2). 

4.2.4.  Correlations  between  the  high  and  low  frequency  parameters 
The  results  of  the  sensitivity  analysis  indicate  a  possible  corre¬ 
lation  between  the  high  and  low  frequency  parameters  at  inter¬ 
mediate  frequencies.  To  clarify  this  point,  the  spectra  in  Fig.  10  were 
fitted  with  a  simple  EEC  that  accounts  for  charge  separation  and 
transport  processes  only  (through  Rhf*  Cdi,c  and  Rc t,c)  which  we  refer 
to  as  the  RRC  equivalent  circuit.  The  parameters  obtained  with  the 
RRC  and  their  corresponding  standard  deviations  are  given  in 
Table  3.  In  order  to  facilitate  the  comparison  with  the  parameters 
estimated  with  the  full  model  including  oxygen  transport,  Table  3 
gives  also  the  relative  difference  AX,  defined  as  |Xrrc  -  XRandies / 
^Randies*  with  X  =  Rhf,  Rc t,c  and  Cdi,o  If  correlations  exist  at  inter¬ 
mediate  frequencies,  the  difference  between  the  estimations  ob¬ 
tained  with  the  two  EEC  could  be  reduced  by  shifting  the  fit  interval 
to  higher  frequencies.  Therefore,  the  parameters  estimations  were 
carried  out  over  three  different  frequency  ranges:  the  maximum 
frequency  is  fixed  at  200  Hz  and  the  minimum  frequency  ^min  is  set 
to  20  Hz,  50  Hz  and  100  Hz.  Indeed,  AX  decreases  for  the  three 
parameters  when  vm[n  is  set  to  50  Hz  instead  of  20  Hz.  However,  the 
following  differences  in  the  values  and  their  variations  can  be 
observed  for  the  three  parameters: 

•  The  high  frequency  resistance  is  only  slightly  affected  by  the 
change  in  the  frequency  range  and  ARhf  remains  below  2%  even 
when  vm[n  is  set  to  20  Hz.  We  can  conclude  that  the  high 
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Table  3 

Mean  cathode  impedance  parameters  and  standard  deviations  estimated  from  a 
series  of  5  impedance  spectra  (Fig.  10)  with  a  simple  EEC  consisting  of  Rh{  in  series 
with  Rctx  and  Cdi,c  in  parallel.  The  table  shows  also  the  normalized  difference  in 
comparison  to  an  estimation  using  the  entire  model  including  oxygen  transport. 


v  [Hz] 

fihf[Q  cm2] 

Rct, c  cm2] 

Cdl,c 

[mF  cm-2] 

AKhf 

[%] 

AKctc 

[%] 

ACdi.c 

[%] 

20-200 

0.156  ±  0.2% 

0.269  ±  0.6% 

22.4  ±0.1% 

1.8 

17 

14 

50-200 

0.151  ±  0.3% 

0.239  ±  1% 

19.9  ±  0.7% 

1.3 

3.9 

0.7 

100-200 

0.148  ±  0.7% 

0.207  ±  3.3% 

18.5  ±  1.6% 

3.2 

9.9 

6.3 

frequency  resistance  can  be  considered  as  independent  of  the 
low  frequency  behavior  of  the  spectra.  This  is  consistent  with 
the  results  of  the  sensitivity  analysis  showing  that  CRhf  is  not 
correlated  with  the  other  sensitivity  coefficients. 

•  ACdi.c  and  ARct,c  depend  strongly  on  the  frequency  range.  The 
impact  of  vmin  is  more  important  on  Rc t)C  than  on  CditC:  ARct,c 
reaches  17%  when  pmin  is  set  to  20  Hz  and  about  4%  when  the 
parameters  estimation  starts  at  50  Hz.  In  comparison,  ACdi>c 
decreases  from  14%  to  0.7%  when  vm[n  is  set  to  50  Hz  instead  of 
20  Hz.  For  an  accurate  estimation  of  Cdi,c  and  Rc t)C  the  use  of  an 
oxygen  transport  impedance  is  thus  necessary. 

Nevertheless,  when  starting  the  parameters  estimation  at  fre¬ 
quencies  above  50  Hz,  AX  remains  below  4%  for  the  three  param¬ 
eters,  confirming  that  the  simple  RRC  equivalent  circuit  is  sufficient 
for  the  estimation  of  the  high  frequency  impedance  parameters  in 
these  conditions.  It  can  be  noted  that  when  vmin  is  further  increased 
to  100  Hz,  AX  increases  again:  this  is  linked  to  the  decrease  of  the 
amount  of  data  available  for  the  parameters  estimation. 

Finally,  it  should  be  mentioned  that  the  results  shown  in  this 
section  are  also  validated  by  sensitivity  analyses  performed  in  other 
operating  conditions  (not  presented  here),  starting  from  experi¬ 
mental  data,  as  well  as  with  series  of  1000  numerically  simulated 
and  noised  spectra. 


5.  Conclusions 

In  this  work,  we  discussed  the  validity  of  the  estimation  of  the 
high  frequency  impedance  parameters  characterizing  the  cathode 
of  a  H2/air  PEMFC,  using  the  Randles  EEC  accounting  for  charge 
separation  and  transport  processes  through  a  high  frequency 
resistance,  a  charge  transfer  resistance  and  a  double-layer  capacity. 
Oxygen  transport  is  accounted  for  by  a  one-dimensional  finite 
Warburg  impedance.  Since  the  mass  transport  limitations  are 
negligible  at  the  anode  when  the  cell  is  fed  with  pure  hydrogen,  it  is 
represented  only  by  a  charge  transfer  resistance  in  parallel  with  a 
capacitance. 

The  parameters  characterizing  the  cathode  reaction  (the  Tafel 
slope  b ,  the  charge  transfer  coefficient  ac  and  the  exchange  current 
density  j0)  can  be  calculated  starting  from  the  charge  transfer 
resistance.  This  allows  us  to  estimate  the  electrode  reaction  pa¬ 
rameters  at  a  given  operating  point  instead  of  giving  only  two  (or 
three)  averaged  values  being  associated  with  the  low  (the  inter¬ 
mediate)  and  the  high  current  density  regimes.  Furthermore,  the 
determination  of  b  and  ac  from  the  Tafel  plot  requires  the  mea¬ 
surement  of  a  whole  current-voltage  curve  and  does  not 
discriminate  between  losses  due  to  mass  transport  and  charge 
transfer.  When  using  EIS  the  interpretation  of  one  single  impedance 
spectrum  is  sufficient  for  the  determination  of  these  parameters, 
which  reduces  drastically  the  measuring  time. 

Complementary  sensitivity  analyses  were  carried  out  in  order  to 
investigate  the  validity  domain  of  the  Randles  EEC  used  for  the 
estimation  of  the  cathode  parameters.  The  numerical  results 


confirm  that  the  cathode  governs  the  cell  impedance:  the  anode 
has  no  significant  impact  on  the  estimation  of  the  cathode 
parameters  and  can  be  neglected.  In  addition,  it  is  possible  to 
determine  the  frequency  range  over  which  the  cathode  parameters 
have  a  significant  impact  on  the  cell  impedance  without  being 
correlated  to  each  other.  Hence,  we  showed  through  the  analysis  of 
the  sensitivity  coefficients  that  Rhf.  Qi.c  and  Rct,c  influence  the 
spectra  rather  in  the  high  and  intermediate  frequency  range, 
whereas  the  diffusion  parameters  of  the  Warburg  impedance,  Rd 
and  id,  have  a  significant  impact  in  the  low  frequency  domain. 
However,  the  results  of  the  sensitivity  analysis  put  forward  the 
existence  of  a  correlation  between  the  charge  transfer  resistance 
and  the  diffusion  parameters  at  intermediate  frequencies.  These 
observations  were  also  confirmed  experimentally  and  the  fre¬ 
quency  range  in  which  the  high  frequency  parameters  can  be 
estimated  independently  of  the  diffusion  impedance  was 
determined. 

Appendix  A.  Cathode  charge  transfer  resistance  at  low 
current  densities 

In  the  sequel,  the  derivation  of  the  cathode  charge  transfer 
resistance  at  low  current  densities  is  presented.  To  start,  the  oxygen 
reduction  reaction  of  the  form  l/202  +  2H+  +  2e-  ->  H20  is 
described  in  a  general  way  by  the  Butler-Volmer  equation  ac¬ 
cording  to 


J  f,c  =  Jo 


ch2o(0)  „„nfac2F 


lh2o 


exp  I  - 


RT 


Vc 


co2(°)exp(z(i_^c 


u02 


(A.1) 


where  ck(x  =  0)  (/<  =  H2 0,02)  stands  for  the  agent  concentrations  at 
the  reaction  interface  at  x  =  0  and  q*  for  their  bulk  concentrations. 
For  very  low  current  densities  (low  overpotentials),  the  exponential 
can  be  approached  by  exp  x  =  1  +  x.  In  addition,  mass  transfer 
limitation  are  negligible  and  the  species  concentrations  at  the  re¬ 
action  sites  can  be  considered  identical  to  their  bulk  values 
ck(0)  =  c£.  The  above  equation  reduces  consequently  to 


(olc2F 


Jf,c  JO  [  RT 


Vc- 


-(1  -  «c)2 F 

RT 


Vc  )  =  ^fVc- 


(A.2) 


Combination  of  (A.2)  with  the  definition  of  the  charge  transfer 
resistance  at  the  cathode  (12)  results  then  in 


Ret 


RT 

Wo ’ 


(A.3) 


showing  that  Rc t  is  independent  from  the  charge  transfer  coefficient 
in  this  current  regime.  This  is  a  general  result  for  RedOx  reactions  in 
electrochemical  systems.  Note  however  that  the  number  of 
exchanged  electrons  might  change  depending  on  the  reaction 
mechanism. 


Nomenclature 


Amacro 

[cm2]  MEA  flat  surface 

Apt 

[cm2]  catalytic  active  area 

MV] 

Tafel  slope 

bio 

[V  dec-1]  Tafel  slope 

Cdi 

[F  cm-2]  double-layer  capacity 

Q 

[mol  cm-3]  concentration  of  species  i 

ci 

[mol  cm-3]  equilibrium  concentration  of  species  i 

Deff 

[m2  s-1]  effective  diffusion  coefficient 
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F0  [V]  open  circuit  voltage  (OCV) 

Eceii  [V]  cell  voltage 

F  [C  mol-1]  Faraday’s  constant,  F  =  96485 
Jceii  [A  cm-2]  cell  current  density 
jo  [A  cm-2]  exchange  current  density 

jc  [A  cm-2]  capacitive  current  density 

jf  [A  cm-2]  faradaic  current  density 

P  [Pa]  gas  pressure 

R  [J  mol-1  K_1]  gas  constant,  R  =  8.314 
Rct  [Q  cm2]  charge  transfer  resistance 

Rd  [0  cm2]  diffusion  resistance 

Fhf  [Q  cm2]  high  frequency  resistance 

Si  [-]  stoichiometric  factor  of  species  i 

T  [K]  temperature 

Z  [Q  cm2] 

impedance 
Greek  letters 

a  [-]  CPE  parameter 

aik  [-]  charge  transfer  coefficient  at  electrode  k 

y  [-]  electrode  roughness 

5  [pm]  diffusion  medium  thickness 

g  [-]  porosity 

%  [V]  activation  overpotential  at  electrode  k 

v  [Hz]  frequency 

id  [s]  characteristic  diffusion  time 

(0  [rad  s-1]  angular  frequency 

Sub-/superscripts 
a  anode 

c  cathode 

Cq2  oxygen  concentration 

W  Warburg 
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